Asymptotics of Universal Probability of Neighboring Level Spacings 

at the Anderson Transition 
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The nearest-neighbor level spacing distribution is numerically investigated by directly diagonalizing 
disordered Anderson Hamiltonians for systems of sizes up to 100x100x100 lattice sites. The scaling 
behavior of the level statistics is examined for large spacings near the delocalization-localization 
transition and the correlation length exponent is found. By using high-precision calculations we 
conjecture a new interpolation of the critical cumulative probability, which has size-independent 
asymptotic form ln/(s) oc —s a with a = 1.0 ± 0.1. 
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The statistical fluctuations in energy spectra of disor- 
dered quantum systems attract at present much attention 
. It is known that by increasing the fluctuations 
of a random potential the one-electron states undergo a 
localization transition, which is the origin of the Ander- 
son metal-insulator transition (MIT) || . The influence of 
the disorder on the wave functions is reflected by the mu- 
tual correlations between the corresponding energy lev- 
els, so that the statistics of energy levels is sensitive to the 
MIT. In the metallic limit the statistics of energy spectra 
can be described by the random-matrix theory (RMT) 
developed by Wigner and Dyson This was shown 

by solving the zero- mode nonlinear <r-model using the su- 
persymmetric formalism ||. Later, perturbative correc- 
tions to the two-level correlation function obtained in the 
RMT were evaluated in the diffusive regime by the im- 
purity diagram technique p0| . In the insulating regime, 
when the degree of disorder W is much larger than the 
critical value W c , the energy levels of the strongly local- 
ized eigenstates fluctuate as independent random vari- 
ables. 

An important quantity for analyzing the spectral fluc- 
tuations is the nearest-neighbor level spacing distribution 
P(s). It contains information about all of the n— level 
correlations. In the metallic regime P(s) is very close to 
the Wigner surmise P\y{ s ) — tts/2 exp (— 7rs 2 /4) jllj (s 
is measured in units of the mean level spacing A). In the 
localized regime the spacings are distributed according 
to the Poisson law, Pp(s) = exp(— s), because the levels 
are completely uncorrelated. The study of the crossover 
of P(s) between the Wigner and the Poissonian limits 
which accompanies the disorder-induced MIT in three- 
dimensional system (3D) was started in Ref. |l^] and 
became the subject of several subsequent investigations 

It was suggested earlier |lj that P(s) exhibits critical 
behavior and should be size-independent at the MIT. In- 
vestigating the finite-size scaling properties of P(s) pro- 
vides not only an alternative method for locating the 



transition 0], but allows also to determine the critical 
behavior of the correlation length [jl4|. A technical ad- 
vantage of the method is that one needs to compute only 
energy spectra and not cigenfunctions and/or the con- 
ductivity. On the other hand, a large number of real- 
izations of the random potential has to be considered. 
In comparison with the well-established transfer-matrix 
method |H| by which one approaches the MIT from the 
localized side, the level-statistics procedure starts from 
the metallic regime. Thus, the two methods can be con- 
sidered to provide complementary information about the 
critical region. 

The suggestion of the existence of a third universal 
level statistics at the MIT excited considerable interest 
in the explicit form of the critical spacing distribution. 
From general considerations for the orthogonal symme- 
try H P(s)(xs at small s. For large s, essentially two dif- 
ferent analytical expressions were proposed |l6[ ]. One of 
them |j],[l2| assumes that P c (s) is a Poissonian for s 3> 1, 
since at the critical point the Thouless energy, which is 
a measure of the number of energy levels that contribute 
to the average conductance of the system, is of order of 
A, while level repulsion is important only for small s. 

A different asymptotic form, P c (s) oc exp(— As a ), was 
proposed (J], by using an analogy between the sequence 
of energy levels and a classical one-dimensional gas of 
interacting fictitious particles. Here a is given by the 
dimensionality d and the localization length exponent v, 



a = 1 + (du) 



(1) 



The result is obtained in the Gibbs model by assuming 
the power law s 2 ~ a for the pairwise interaction between 
the particles ||. The latter distribution decays faster 
than the Poissonian (a = 1), but slower than the Wigner 
surmise (a = 2). Several numerical calculations for the 
3D Anderson model were recently performed PJT^| in or- 
der to analyze P c (s). The results were found to be consis- 
tent with the latter of the above suggestions with an ex- 
ponent a w 1.2 — 1.3 [y « 1.5). However, since the round- 
ing errors in the calculations for large s are such that 
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TABLE I. Numerical parameters for various cube sizes L 
at W C =1QA. M: number of samples, N s : total number of 
spacings, A: mean level spacing, p = (AL 3 )" 1 : density of 
states; a and A c : quantities of Eq. (jij). All levels lie within 
the energy interval \E\ < 4.45. 
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a = 1 cannot be completely ruled out, the asymptotic 
form of P c (s) is still an open question, and the subject of 
presently on going and controversial discussions. In this 
Letter we present the results of detailed high-precision 
numerical investigations of the critical level spacing dis- 
tribution. Our findings solve the above controversy. Pre- 
liminary results have been published previously fl8[| . 

By diagonalizing the Anderson Hamiltonian with a 
Lanczos algorithm Jl4| specifically modified for systems 
containing up to 10 6 lattice sites, which were not achieved 
in previous works, we examined both the critical be- 
havior and the finite-size scaling properties of the inte- 
grated probability distribution of neighboring spacings 
I(s). Our main result is that the asymptotic form of 
critical I c (s) and, therefore, P c {s) at large s is very close 
to a Poissonian decay, as the leading term, thus confirm- 
ing the ideas of @JT^]. In addition, by using the size- 
independence of I c (s) at the MIT and investigating the 
scaling of I(s) with the system size L and W, we estimate 
the correlation length exponent v. 

The Anderson model |l9| is defined by H = 

E„ £ n4 a n + E„/ ro (4 a m + C - C 0, wher e 4 (flj is the 

creation (annihilation) operator of an electron at a site 
n, with m denoting the nearest neighbors of n. The site 
energies e n are measured in units of the overlap integral 
between adjacent sites. They are independent random 
variables that are distributed around e = according to 
a box distribution of width W. A simple cubic lattice 
with periodic boundary conditions was used. We com- 
puted the electron spectra of cubes of linear size ranging 
from L = 5 to 100 for various W. It is known from 
the transfer-matrix method pfj, that in the center of 
the band (e = 0) W c « 16.4. The spectrum was prop- 
erly "unfolded" by fitting the integrated density of states 
around e = to polynomial splines. The numerical re- 
sults at the MIT are summarized in Table |. It should be 
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FIG. 1. Level spacing distribution P(s) for various sys- 
tem sizes at the critical disorder W c . Dashed-dotted line is 
Pakl(s). Full line is derivative of I c (s) from interpolation 
formula (Q). 

noted that the numerical diagonalization of giant sparse 
matrices of order of 10 5 -10 6 is highly nontrivial. 

Fig. shows P(s) calculated at the MIT. As expected, 
it is L-independent. To cover the whole range of spacings, 
the interpolation formula Pakl{s) — Bs exp(-is a ) 
has been proposed in Due to normalization A = 
[r(3/cv)/r(2/a)] Q and B = aA 2 / a /T{2/a). The best fit 
using the x 2_ criterion in the interval < s < 4 yields 
a = 1.48 ± 0.08 with a confidence level 0.95. The fitted 
exponent a is markedly larger than that given by (Q). 
For s > 3 one observes an increasing deviation between 
Pakl(s) and the computed histogram. This shows that 
fitting near s ~ 1 does not provide reliable information 
about a, because the exponential tail of P(s) contributes 
to the relative accuracy only with a very small weight. 
Therefore it is imperative to investigate the asymptotic 
behavior at large s, not including data from the region 
< s < 2. 

In what follows, wc consider the cumulative level spac- 
ing distribution function I(s) = P(s')ds'. It gives 
the probability to find neighboring energy levels with a 
separation E > s A. The integration does not change 
the asymptotic exponential behavior of P(s). Since 
s > 0, 7(0)=1, and by normalization to the total num- 
ber of spacings in a given interval, / °° I(s) ds = 1. The 
Wigner surmise and the Poisson distribution yield 
Iw{s) = exp(— 7rs 2 /4) and Ip{s) — exp(— s), respec- 
tively. The numerical evaluation of I(s) is similar to that 
of the density of states in unfolding the spectrum. By 
arranging the spacings in a descending sequence one can 
very accurately construct the histograms of I(s) 

Using the common statistical hypothesis at large s 

\nI c (s) = -A c s a , (2) 
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FIG. 2. Critical probability of neighboring spacings I(s). 
Solid line is Eq. (j^) with a = 1. Dashed-dotted line is power 
law with q = 1.24 from Eq. (|l|). Dashed and dotted line are 
Iw(s) and Ip(s), respectively. 

we calculated A c and a for various L (see Table |) . Inde- 
pendently of L the result is a = 1.0±0.1. The numerical 
data of In I c (s) shown in Fig. || are better described by 
a linear law for s > 3, so that I c {s) oc exp(— ^4 c s) with 
A c = 1.9 ± 0.1. This is similar to the insulating regime, 
although the decay rate A c is larger than unity due to 
the level repulsion. The power law with the exponent 
a w 1.2, which was recently obtained jljj ] by an analysis 
of the shape of P(s) in the range < s < 5 for system 
sizes L < 21 deviates from our results for s > 4. 

The linear asymptotic behavior of ln/ c (s) is in con- 
trast to the power law with a w 1.31 obtained numeri- 
cally H for smaller systems L < 12. The reason for this 
discrepancy is the following. The energy interval E con- 
sidered in H is so narrow that it contains only ten spac- 
ings on the average, that results in a cut-off of P c {s) at 
s f» 10 A. Thus, some fraction of the spacings s < 10A is 
not taken into account, causing the faster decay of P c (s). 
In our calculations the interval is wide enough, covering 
approximately half of all of the eigenvalues. However, 
such a choice of E does not lead to the undesirable mix- 
ture of the extended and the localized states. This is due 
to a peculiarity of the box distribution of the site ener- 
gies e n . It follows from the localization phase diagram 
{W C ,E C } [ ^5|j22| , that the critical disorder W c is almost 
independent of the energy when \E C \ < 6. In order to in- 
vestigate how the width of the energy interval influences 
the level statistics, we calculated J(s) for E/A = 10 2 , 10 3 
and 10 4 , provided that all levels satisfy the critical con- 
dition L < £(e) oc \e/E c — l\~ v . The results were practi- 
cally the same within the statistical uncertainties. This 




FIG. 3. Probability I(s) for L=6 and 28 at W=12, 14, 15, 
16, 16.4, 17, 18 and 20 shown consecutively from the left to 
the right. Dashed (dotted) line is the Wigner (Poisson) limit. 

implies the equivalence of averaging over the spectrum 
and over the random potential. Indeed, due to dimin- 
ishing the spacing with the size A oc L 3 , the averaging 
for smaller cubes is performed over many samples, while 
for our largest systems L=80 and 100 a single realiza- 
tion without ensemble averaging is even sufficient to get 
similar distributions with comparable precision. For nu- 
merically describing a crossover between small and large 
s, we propose an explicit form of the new interpolation 
function 

I c (s) = exp[ M - vV + (A c s) 2 } (3) 

with a coefficient fi « 2.21. Although we do not provide a 
rigorous analytical prove, it gives the excellent fit all over 
the range of the computed spacings. The corresponding 
P(s) shown in Fig. [I] fulfills the both normalization con- 
ditions. 

To study the finite-size scaling behavior of I(s) for 
large s, we extended the calculations to other degrees 
of the disorder W close to W c for various system sizes. 
The calculations were performed for an ensemble of dif- 
ferent samples. The number of samples for each given 
pair of L and W was chosen such that N s ~ 10 5 spac- 
ings were obtained. We have also carefully checked the 
sensitivity of the results to the number of realizations. 
No change was observed within the error bars when in- 
creasing the system size on the expense of the number of 
realizations and vice versa. By increasing W the spacing 
distribution for fixed L changes continuously from Iw(s) 
to Ip(s) (Fig. ||). The steepness of the crossover depends 
on L. For larger sizes I(s) changes faster between the two 
limiting regimes. At W c ~ 16.4 the spacing distribution 
has almost the same asymptotic form for all L from 5 to 
100. This reflects the universality of the level statistics 
exactly at the MIT 0. 

For finite L the distribution I(s) exhibits scaling in 
the vicinity of W c . Within the critical region, L < £(W), 
it is reasonable to assume that the linear slope of ln/(s) 
is governed by the one-parameter scaling law, A(W, L) = 
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FIG. 4. Scaling variable A as a function of the disorder W 
for different L, showing critical behavior near the MIT. Inset: 
the one-parameter dependence of A on L/£(W). 

f(L/£(W)). Fig. H shows the disorder dependence of A 
near the critical point for various L. A(W C ) does not de- 
pend on L. By introducing a scaling parameter, the cor- 
relation length £(W) @H, we found a common scaling 
curve consisting of two branches corresponding to the 
delocalized and the localized regimes for A > A c and 
A < A c , respectively, as shown in the inset of Fig. ||. 
The critical exponent v was determined in a similar way 
as previously where only the small- s part of P(s) 
was used. We found v = 1.4 ± 0.15 in agreement with 
the result obtained earlier by completely different meth- 
ods|,[5@. 

In conclusion, we present the first large-scale numer- 
ical results on the statistics of the energy levels near 
the disorder-induced MIT for systems of sizes up to 
L 3 = 100 3 sites. A comparative analysis with results 
obtained from various analytical approaches and other 
numerical studies is performed. At the critical point the 
asymptotic universal probability of energy level spacings 
has a Poisson-like form I c (s) oc exp(— A c s). We believe 
that the simple exponential asymptotics of the critical 
level spacing distributions are valid not only for the or- 
thogonal symmetry (with spinless electrons and without 
magnetic field), but also for other universality classes: 
the unitary (in the presence of the magnetic field) and 
the symplectic (in the presence of spin-orbit coupling) 
classes. Recent computer simulations |23] corroborate 
that the decay rate A c is almost insensitive to the fun- 
damental symmetry. However, it could depend on the 
physical dimensionality. Finally, we have determined the 
influence of the disorder of the system on the exponen- 
tial tail of I(s), and constructed numerically the corre- 



sponding scaling function. The critical exponent of the 
correlation length was calculated, v « 1.4. 
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